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Abstract 

A detailed study is made of the effects of vari- 
ations in lamination and material parameters of 
thin-walled composite frames on their vibrational 
characteristics. The structures considered are semi- 
circular thin-walled frames with I and J sections. 
The flanges and webs of the frames are modeled 
by using two-dimensional shell and plate finite ele- 
ments. A mixed formulation is used with the fun- 
damental unknowns consisting of both the general- 
ized displacements and stress resultants in the frame. 
The frequencies and modes predicted by the two- 
dimensional finite-element model are compared with 
those obtained from experiments, as well as with the 
predictions of a one-dimensional, thin- walled-beam, 
finite-element model. A detailed study is made of 
the sensitivity of the vibrational response to varia- 
tions in the fiber orientation, material properties of 
the individual layers, and boundary conditions. 

Introduction 

The physical understanding and the numerical 
simulation of the dynamic response of laminated 
anisotropic structures have recently become the fo- 
cus of intense efforts because of the expanded use 
of fibrous composites in the aerospace, automotive, 
shipbuilding, and other industries, and because of 
the need to establish the practical limits of the dy- 
namic load-carrying capability of structures made 
from these materials. Experimental studies have 
been conducted on the free vibration and impact re- 
sponse of thin- walled composite frames and stiffeners 
(e.g., see Boitnott et al. 1987; Boitnott and Fasanella 
1989; Collins and Johnson 1989; and Chandra, Ngo, 
and . Chopra 1988). One-dimensional theories have 
been developed for the static, vibration, and buckling 
analyses of thin-walled-frame structures (e.g., Vlasov 
1961; Gjelsvik 1981; Nowinski 1966; and Panovko 
and Beilin 1969). However, no systematic assessment 
has been made of the range of validity of the basic 
assumptions of these theories. Approximate analyt- 
ical and numerical techniques have been applied to 
the study of the vibrational response of isotropic and 
composite stiffeners (e.g., see Hasan and Barr 1974; 
Vermisyan and Galin 1972; Rao 1975; Vasilenko and 
Trivailo 1980; Narayanan, Verma, and Mallik 1981; 
Ali 1984; Gupta, Venkatesh, and Rao 1985; Potiron 
et al. 1985; Riickschloss 1985; Wekezer 1987; Reh- 
field, Atilgan, and Hodges 1990; Stemple and Lee 
1988; and Bishop, Cannon, and Miao 1989). Few 
publications exist in which the effects of variations 
in lamination and geometric parameters of composite 
panels on their vibrational characteristics are studied 
(see Teh and Huang 1980 and Bank and Kao 1989). 


However, none of these publications consider thin- 
walled composite frames. 

The present study is an attempt to fill this void. 
Specifically, the objective of this paper is to sum- 
marize the results of a recent study on the effects 
of variations in the lamination and of geometric pa- 
rameters of thin-walled composite frames on their 
vibrational characteristics (frequencies, and energy 
components associated with different modes). The 
frames considered are semicircular, are made of thin- 
walled graphite-epoxy material with I and J sections, 
and have a 36-in. radius (see fig. 1). 


Symbols 


A 

cross-sectional area of one- 
dimensional- beam model 

Ay, A z 

effective shear areas for one- 
dimensional-beam model in y- 
and z-directions, respectively 

A n 

extensional stiffness of laminate 
(flanges or web) in xi-direction 

^33 

in-plane shear stiffness of 
laminate 

Bu, 

bimoment in one-dimensional- 
beam model 

b, h,b2, 

& 3>&4 

flange (or skin) dimensions (see 
tables 2 and 3) 

Ci 

multipliers (see eqs. (5) and (7)) 

d, d\ , d2 

web dimensions (see tables 2 
and 3) 

E,G 

effective Young’s and shear 
moduli of equivalent isotropic 
material 

E l ,E t 

elastic moduli of individual layers 
of laminate (flanges or web) in 
direction of fibers and normal to 
it, respectively 

[F] 

matrix of linear flexibility coeffi- 
cients for an individual element 

Glt, Gtt 

shear moduli in plane of fibers 
and normal to it, respectively 

w 

vector of stress resultant (or 
internal force) parameters 

h 

total thickness of laminate 

hi,h.2, h$, 
b-A% h{), 

h 7 

wall thicknesses (see tables 2 
and 3) 



lyi Iz > lyz 

lu 

J 

K 

m 

i 

M y ,M z ,M t 

Mi2,M 2 i 

[MUM]; 

N2,N\2- 

NL 

N x 

TOQ] * 

Qy j Qz 
Qi,Q2 
R 

Ri 

[5] 

u c 

U t f,U Wi U bf 


second moments of cross section 
(moments and product of inertia) 
of one-dimensional-beam model 

principal second sectorial mo- 
ment of cross section of one- 
dimensional-beam model 

Saint- Venant torsion constant of 
cross section of one-dimensional- 
beam model 

kinetic energy 

generalized stiffness matrix for an 
individual element (see eqs. (3)) 

length of individual finite element 

bending and twisting moments in 
one-dimensional-beam model 

bending stress resultants in two- 
dimensional model 

consistent and generalized 
mass matrices for an individual 
element (see eqs. (3)) 

extensional stress resultants in 
two-dimensional model 

total number of layers 

axial force in one-dimensional- 
beam model 

matrices associated with con- 
straint condition and regulariza- 
tion term in the functional for 
one-dimensional-beam model 

transverse shear forces in one- 
dimensional-beam model 

transverse shear stress resultants 
in two-dimensional model 

radius of curvature of center- 
line of frame (used in one- 
dimensional-beam model) 

outer radius of curvature of frame 
(see fig. 1) 

strain-displacement matrix for 
individual element 

total complementary strain 
energy of frame 

contributions of top flange, web, 
and bottom flange (including 
skin) to total complementary 
strain energy 


Ui,u 2 

complementary strain-energy 
components associated with in- 
plane and out-of-plane forces, 
respectively 

U 3 

complementary strain-energy 
component associated with forces 
neglected in one-dimensional- 
beam model 

U,V,W 

displacement components in 
coordinate directions for one- 
dimensional-beam model 


axial and transverse displace- 
ments of one-dimensional-beam 
model at y = z = 0 

U\,U 2 ,W 

displacement components of two- 
dimensional model in xi,x 2 ,X 3 
coordinate directions 

u'i,u' 2 ,w' 

displacement components of two- 
dimensional model in x'^x^Xg 
coordinate directions 

w 

vector of nodal displacements 

x,y,z 

centroidal orthogonal coordinate 
system used for one-dimensional- 
beam model 

Xl,X2,X 3 

local orthogonal coordinate 
system used in conjunction with 
two-dimensional model (for the 
web and each of the two flanges) 

®2> ®3 

global Cartesian coordinates used 
for two-dimensional model 

m 

vector of element degrees of 
freedom 

{Z} 

particular solution (see eqs. (5) 
and (7)) 

Ixyilxzi 
7:n/5 7 xz 

transverse shear strains in 
one-dimensional-beam model (see 
eqs. (A2» 

e 

penalty parameter 

£x 

extensional strain in one- 
dimensional-beam model 

£ 0 

extensional strain of centerline of 
one-dimensional-beam model 

6 

angle that a typical cross section 
of frame makes with x [ , x 2 plane 
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0° rate of twist of one-dimensional- 

beam model 

Ky, ft®, k® curvature changes and twist of 
one-dimensional-beam model 

{A} vector of Lagrange multiplier 

parameters 

A Lagrange multiplier 

A* lamination and material 

parameters 

i/LT major Poisson’s ratio of 

individual layers 

o X: &xy, &xz normal and shearing stresses on 
cross section of beam 

7r, n HR functionals defined in equa- 

tions (A7) and (A8) 

p mass density of material 

<£l, 4>2 rotation components of two- 

dimensional model referred to 
local coordinate system X\,X 2 

(j>l,<f> 2, <t> 3 rotation components of two- 

dimensional model referred 
to global coordinate system 
x' v x' 2 ,x' 3 

(f>y, <fP z rotation components in one- 

dimensional-beam model 

\P° strain parameter in one- 

dimensional-beam model 

lj frequency of vibration 

u sectorial coordinate (warping of 

cross section for a unit rate of 
twist) 

d = d/dx 

ID one-dimensional-beam model 

2D two-dimensional model 

Subscript: 

s shear center 

Superscript: 

t matrix transposition 

Analysis 

Computational Models 

Two computational models are used for the thin- 

walled composite frames considered in the present 


study. In the first model, the flanges and web are 
modeled by using two-dimensional shell and plate 
finite elements. The second model is a finite-element 
discretization of the one-dimensional Vlasov type 
thin-walled-beam theory. Herein, the two models 
are referred to as two-dimensional (2D) and one- 
dimensional (ID) finite-element models, respectively. 

Mathematical Formulation 

Two-dimensional models . The analytical for- 
mulation for the two-dimensional models is based 
on the Sanders-Budiansky shell theory with the ef- 
fects of transverse shear deformation, and laminated 
anisotropic material response included. A mixed 
formulation is used in which the fundamental un- 
knowns consist of the generalized displacements and 
the stress resultants in the frame. (See fig. 2 for the 
sign convention.) 

Bicubic shape functions are used to approximate 
each of the generalized displacements and stress re- 
sultants. There are 16 displacement nodes and 128 
stress-resultant parameters in each element. The 
stress resultants are allowed to be discontinuous at 
interelement boundaries. The element characteristic 
arrays are obtained by using the two-field, Hellinger- 
Reissner, mixed-variational principle. 

One- dimensional models . The analytical for- 
mulation for one-dimensional models is based on a 
form of the Vlasov thin- walled-beam theory with the 
effects of flexural-torsional coupling, transverse shear 
deformation, and rotary inertia included. The fun- 
damental unknowns consist of seven internal forces 
and seven generalized displacements of the beam (see 
fig. 3 for the sign convention). The element char- 
acteristic arrays are obtained by using a modified 
form of the Hellinger-Reissner mixed variational prin- 
ciple. The modification consists of augmenting the 
functional of that principle by two terms: (1) the 
Lagrange multiplier associated with the constraint 
condition relating the rotation of the cross section 
and the twist degrees of freedom, and (2) a reg- 
ularization term that is quadratic in the Lagrange 
multiplier. Only C° continuity is required for the 
generalized displacements. Lagrangian interpolation 
functions are used for approximating each of the gen- 
eralized displacements, internal forces, and Lagrange 
multiplier. The polynomial functions for the inter- 
nal forces and the Lagrange multiplier are one degree 
lower than those of the generalized displacements. In 
the present study, quadratic polynomials are used in 
approximating the generalized displacements. Linear 
polynomials are used in approximating each of the 
internal forces and the Lagrange multiplier. The in- 
ternal forces and the Lagrange multiplier are allowed 
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to be discontinuous at interelement boundaries. For 
each element, the total number of generalized dis- 
placement parameters is 21, the total number of in- 
ternal force parameters is 14, and the total number of 
Lagrange multiplier parameters is 2. The fundamen- 
tal equations of the thin-walled-beam theory used in 
the present study are given in Noor, Peters, and Min 
(1989) and are summarized in the appendix. 

For quasi-isotropic laminated composites, numer- 
ical experiments to be described subsequently have 
demonstrated that reasonably accurate results can 
be obtained with the one-dimensional model when 
the laminated composite is replaced by an equiva- 
lent isotropic material with the following Young’s and 
shear moduli: 


Array 


Two-dimensional models 


One- dimensional models 


{Z} 




[*] 



r-F s 
s 1 

_ ■ Q l 


Q 




[ifr] 




E = An/h 
G = A&/h 


(1) Sensitivity of Vibrational Response to 

(2) Variations in Lamination and Material 
Parameters 


where An and A33 are the extensional stiffness in 
the ^-direction and the in-plane shear stiffness used 
in the classical lamination theory, and where h is 
the total wall thickness (of the flange or web). This 
approximation was adopted in the present study. 


Finite-Element Equations 

The finite-element equations for each individual 
element of the two-dimensional and one-dimensional 
models can be cast in the following compact form: 

([k]-u?[M]){Z} = 0 (3) 


where {Z} is the vector of the element degrees of free- 
dom, u) is the frequency of vibration, and [k] and [ M ] 
are the generalized stiffness and mass matrices. The 
forms of {Z}, [AT], and [m] are defined in the follow- 
ing table, where {H }, {Y}, and {A} are the vectors of 
stress-resultant (or internal force) parameters, nodal 
displacements, and Lagrange multiplier parameters, 
respectively; [F] is the matrix of linear flexibility co- 
efficients; [S] is the strain displacement matrix; [P] 
and [Q] are matrices associated with the constraint 
condition and the regularization term in the func- 
tional, respectively; [M] is the consistent mass ma- 
trix; e is a penalty parameter associated with the 
regularization term; superscript t denotes transposi- 
tion; and a dot (•) refers to a zero submatrix. The 
explicit form of the arrays in the following table is 
given in Noor and Andersen (1982) and Noor and 
Peters (1983) for the two-dimensional models and in 
Noor, Peters, and Min (1989) for the one-dimensional 
models: 


The expressions for the sensitivity derivatives of 
the frequency and response vectors with respect to 
the lamination and material parameters A^ of the 
composite frames are given by (Nelson 1976) 


8U 2 _ yy 

d\ i 

Elements 




d[M\ 

d\, 


{ z } (4) 


and 


^ = {Z} + Ci{Z} (5) 

where {z} represents a particular solution of the 
equations 


([*]-“ 2 i")>^r= 

(d[K] 2 d[M] aw 2 , . A 

1«r" W'W' I J ,Z1 (6) 

and Ci are multipliers given by 

« = - E (wimzy+izy'l^fliz}) 

Elements V l U 

( 7 ) 

In equations (4) to (7), the eigenvectors are as- 
sumed to be normalized with respect to [M ]; that 
is, 

{ZY[M]{Z} - 1 (8) 

The expressions for the total complementary 
strain energy of the frame U c and its derivatives with 
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respect to A* axe given by 


and 


u c = \ E (9) 

Elements 



For the purpose of obtaining analytic derivatives with 
respect to some of the lamination parameters, such 
as the fiber orientation angle of different layers, it 

is convenient to express in terms of as 

follows: 


d[F] 

d\i 


= -[F] 


m - 1 

d\i 




( 11 ) 


The matrix is evaluated using the analytical 

derivatives of the material stiffness matrix of each 
laminate (flanges and web). The material stiffness 
matrix of the laminate is given in Jones (1975). 


Experimental and Numerical Studies 

Apparatus and Test Procedure 

Specimens. Two specimens, an I-section and a 
J-section frame (fig. 4), were tested in the present 
study. Nominal dimensions of each cross section are 
given in figure 1. Weights of the frame sections were 
3.181 lb and 4.085 lb for the I and J frames, respec- 
tively. The frame sections were made from AS4/5208 
graphite-epoxy unidirectional tape laid up in a man- 
ner that resulted in essentially uniform stiffness prop- 
erties in the circumferential direction (i.e., the stiff- 
ness coefficients are independent of 6). The ma- 
terial properties for the individual layers are given 
in figure 1. The laminate stacking sequence was 
[±45/0/90] s for the I-section and [±45/0/90] 2s for 
the J-section. Each frame section was semicircu- 
lar with a diameter of 72 in. Bonded to the out- 
side flange of each frame was a 16-ply [±45/0/90]2 5 , 
quasi- isotropic skin made of the same material. The 
frame sections were constructed so that the skin 
would extend 0.5 in. beyond each side of the bottom 
flange of the frame. The nominal dimensions of the 
I- and J-section frames are given in figure 1, and the 
actual (measured) dimensions are given in tables 1 
and 2. 

Apparatus and procedure . Figure 5 is a 

schematic of the experimental setup, and figure 4 is a 
photograph of the setup and specimens. The ends of 


the frame sections were potted in a fixture that was 
bolted to a large steel-beam backstop. (See fig. 4.) 

An air shaker, connected to an air compressor, 
was used to excite all test specimens. Excitation was 
both in plane (radially) and out of plane. For in- 
plane excitation, the shaker was positioned so that 
the pulses of air struck approximately normal to the 
surface of the skin. For out-of-plane excitation, a 
piece of Dow Chemical Co. Styrofoam was attached 
to the side of the frame by double-sided adhesive 
tape. Pulses of air struck the flat face of the Sty- 
rofoam normal to the face. The position of the air 
shaker was adjusted when the excitation was striking 
on a node. 

A miniature accelerometer was attached at a fixed 
location to the frame sections with double-sided ad- 
hesive tape. Output from the accelerometer was am- 
plified and displayed along the vertical axis of an os- 
cilloscope. Natural modes were determined by tuning 
the excitation frequency of the air shaker to produce 
a maximum acceleration of the vertical deflection on 
the oscilloscope. Output also passed through a low- 
pass filter and was displayed as vibrational frequency 
on a frequency counter. 

A hand-held velocity probe was moved along the 
frame to determine node locations and mode shapes. 
The output of the probe was displayed along the 
horizontal axis of the oscilloscope. The probe and 
accelerometer outputs combined to create a Lissajous 
pattern on the oscilloscope. A phase shift in the 
Lissajous pattern occurred when the velocity probe 
passed over a node. 

The nodal locations were mapped manually dur- 
ing the vibration survey of the frames. Consequently, 
the only nodal lines monitored were those associ- 
ated with gross in-plane and out-of- plane motions. 
Other nodal lines, associated with localized defor- 
mation patterns, were not surveyed. These localized 
deformations were noticeable in some of the higher 
vibration modes, with complex deformation patterns 
and/or strong coupling between in-plane and out-of- 
plane motions. 

Finite-Element Grids 

Two-dimensional models were generated for the 
actual frames (test specimens) described in the pre- 
ceding section and for the corresponding frames with 
nominal dimensions. Herein, the frames with actual 
and nominal dimensions will be referred to as the 
actual and nominal frames, respectively. For the ac- 
tual frames, spline interpolations were used to gener- 
ate the wall thicknesses and coordinates of the nodal 
points. Isoparametric finite elements were used to 
approximate the variations in stiffness and geome- 
try. The one-dimensional models considered herein 
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are for the frames with nominal dimensions. The 
grids used for both the one-dimensional and two- 
dimensional models are described subsequently. 

Two-dimensional models . An 18 x 8 grid was 
used for modeling the whole I-section frame. In this 
grid, two elements were used to model each of the 
web, top flange, and bottom flange sections. The 
part of the skin adjacent to the bottom flange section 
was treated as part of the flange. One element was 
used to model each of the two parts of the skin section 
that extended beyond the bottom flange. (See fig. 1.) 
The middle surfaces of the top flange and the web 
were taken to be their reference surfaces. The middle 
surface of the combined bottom flange and skin was 
taken to be the reference surface. 

An 18 x 7 grid was used for modeling the whole 
J-section frame. The distribution of the elements was 
similar to that for the I-section frame. Only one 
element was used to model the top flange section. 
(See fig. 1.) 

Totally clamped and partially clamped support 
conditions were considered. For totally clamped 
supports, all six generalized displacements were re- 
strained (u[ = u f 2 = w f = 0^ = <f > 2 = <£3 = 0). The 
partially clamped conditions were obtained from the 
totally clamped case by successively removing the re- 
straints on one, as well as on combinations, of the 
displacement and rotation components. 

One- dimensional models . A uniform grid 

of 24 elements was used in modeling each of the 
I-section and J-section frames. The principal sec- 
torial properties of the cross section were evaluated 
with the Fortran program listed in Coyette 1987. 

Identification of Modes and Estimation of 

Error in One-Dimensional-Model 

Predictions 

The two-dimensional models can be used to iden- 
tify the in-plane, out-of-plane, and coupled modes 
and to estimate the error in the predictions of the 
one-dimensional models. These objectives can be 
accomplished by decomposing the complementary 
strain energy U c (eq. (9)) associated with each vi- 
bration mode into three components, C/2 , and 
C/3 (see table 3). The first two components, U\ 
and C/2, are associated with the in- plane and out- 
of-plane stress resultants, respectively. The third 
component, C/3, is associated with the stress resul- 
tants that are peculiar to two-dimensional plates and 
shells (not present in one-dimensional-beam models). 
The in-plane and out-of-plane modes correspond to 
the modes for which U\/U c and U 2 /U c are close to 
1, respectively. The strongly coupled modes corre- 
spond to nearly equal values of U\/U c and U 2 /U c . 


The ratio Us/U c is indicative of the error in the one- 
dimensional-model predictions. 

It is also useful to partition the total complemen- 
tary strain energy associated with each mode into 
three components, U t f : U w , and C/&y; these compo- 
nents represent the contributions of the top flange, 
web, and bottom flange (including the skin). 

Comparison of Experimental and 

Finite-Element Results 

The results of the experimental and numeri- 
cal studies are summarized in figures 6 to 10 and 
table 4 for the I-section frame, and in figures 11 to 15 
and table 5 for the J-section frame. Figures 6(a) and 
11(a) are bar charts for the experimental frequencies 
and the frequencies obtained by the two-dimensional 
finite-element model for the actual I-section and 
J-section frames, respectively. For the finite-element 
model, three cases are considered — totally clamped 
edges (with both translational and rotational re- 
straints), partially clamped edges with 0 f 2 not re- 
strained, and partially clamped edges with u[ in the 
flanges and 02 not restrained. 

The maximum and minimum values of the fre- 
quencies obtained by the two-dimensional finite- 
element model (corresponding to the totally clamped 
and partially clamped edges) are shown in fig- 
ures 6(b) and 11(b), along with the experimental fre- 
quencies. (See also tables 4 and 5.) The experimen- 
tal frequencies associated with modes 9 and 10 of the 
I-section frame, and with modes 9, 10, and 11 of the 
J-section frame, respectively, are close in frequency 
and have very close nodal locations. Henceforth these 
modes will be referred to collectively as mode 9. Also, 
the 12th mode of the I-section frame (table 4) was 
missed in the experimental survey, which is indicative 
of the difficulty of determining the high-frequency 
modes. The nodal locations of the succeeding ex- 
perimental frequency for both the I- and J-section 
frames are close to those of the finite-element model. 
The fact that only one of the multiple experimental 
frequencies with close nodal locations is predicted by 
the finite-element model may be attributed to imper- 
fections in lamination and material properties and/or 
to geometric nonlinearities that were not incorpo- 
rated into the finite-element model. Figures 6(c) 
and 11(c) are bar charts for the frequencies obtained 
by two-dimensional models of the actual and nomi- 
nal frames along with those of the one-dimensional 
model. 

Figures 7 and 12 are bar charts of the two decom- 
positions of the complementary strain energies, asso- 
ciated with the different vibration modes, described 
in the preceding section. The ordinates in figures 7(a) 
and 12(a) represent the ratios of U\/U c , U 2 /U c , and 
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Us/U c , and the ordinates in figures 7(b) and 12(b) 
represent the ratios of U t f/U c , U w /U c , and U^j/U 0 
for each of the modes. 

The mode shapes associated with the first 10 fre- 
quencies are shown in figures 8 and 13. Three views 
are shown for the deformations associated with each 
mode — side view, top view, and end view. Also 
shown are the nodal lines of the w 9 displacement 
on the top and bottom flanges. As can be seen 
in figures 8 and 13, the deformation patterns asso- 
ciated with higher modes are fairly complex. As 
mentioned previously, the only experimental nodal 
lines monitored are those associated with gross in- 
plane and gross out-of-plane motions. Generally, 
good agreement between the finite-element and ex- 
perimental nodal lines is observed in these cases. 
Other nodal lines, associated with localized deforma- 
tions, are shown only for the finite-element solutions. 

The sensitivities of the vibration frequencies to 
the fiber orientation angles of the top flange, web, 
and bottom flange and skin are depicted in figures 9 
and 14. The ordinates in figures 9 and 14 represent 
the sensitivity derivatives with respect to the indi- 
cated fiber angles. Each of the sensitivity derivatives 
is normalized by dividing it by the corresponding fre- 
quency of vibration. The sensitivities of the vibra- 
tion frequencies to the material parameters E £, 

Git > an d Gtt are shown in figures 10 and 15. The 
ordinates in figures 10 and 15 represent the sensi- 
tivity derivatives with respect to the indicated elas- 
tic moduli. Each of the sensitivity derivatives is di- 
vided by the corresponding frequency and multiplied 
by the corresponding elastic modulus. The effects of 
boundary conditions on the frequencies obtained by 
the two-dimensional finite-element model are shown 
in tables 4 and 5. 

An examination of the experimental and finite- 
element results (figs. 6 to 15 and tables 4 and 5) 
reveals the following: 

1. Reasonably good correlation is observed be- 
tween numerical simulation and experiment for the 
I-section frame (fig. 6(a)). The ratios of the first 
five experimental frequencies to the corresponding 
finite-element frequencies ranged from 0.90 to 1.00 
(table 4). For the J-section frame, the correlation is 
not as good (fig. 11(a)). The corresponding ratios 
for the first five frequencies were from 0.86 to 1.01 
(table 5). 

2. Most of the experimental frequencies for the 
I-section frame and the J-section frame are between 
those for the totally and partially clamped supports 
(with both u[ in the flanges and (p 2 not restrained), 
especially for the higher modes. For some of the 
modes, the experimental frequencies are closer to 
the partially clamped support case (e.g., modes 5, 


7, 10, 12, and 13 (fig. 11(b))). For the I- and 
J-section frames, the finite-element model predicted 
only one of the multiple experimental modes with 
close nodal lines. The other experimental frequencies 
were between those for the totally and partially 
clamped supports (with both in the flanges and 
<f > 2 not restrained (figs. 6(b) and 11(b))). 

3. The lowest five frequencies obtained by the 
one-dimensional model are reasonably close to those 
obtained by the corresponding two-dimensional 
model, especially for the J-beam, where the errors 
in the predictions of the one-dimensional model were 
well below 10 percent. (See figs. 6(c) and 11(c).) 

4. Identification of the modes as in plane or out 
of plane can best be accomplished by examining the 
energy components, U\/U c and U^/U 0 , associated 
with the in-plane and out-of-plane forces, respec- 
tively (figs. 7(a) and 12(a)). Also, the minimum er- 
ror to be expected when using one-dimensional thin- 
walled beams can be estimated by computing the 
ratio of the energy associated with the forces ne- 
glected in thin-walled beams to the total energy 
U%/U c . (See figs. 7(a) and 12(a).) 

5. The coupling between in-plane and out- 
of-plane deformations is more pronounced in the 
J-section frame than in the I-section frame. For ex- 
ample, the first 20 modes for the I-section frame had 
either U\/U c or U<ilU c > 0.75. On the other hand, 
only modes 1 to 4, 6, 8, and 10 in the J-section 
frame had U\/U c or U^/U 0 > 0.75. For the higher 
modes, neither U\/U c nor C/ 2 /U c was close to 1. (See 
figs. 7(a) and 12(a).) 

6. For the I-section frame, the contributions of 
the top and bottom flanges to the total energy asso- 
ciated with different modes far exceeded that of the 
web. The ratio of the strain energy in the web to the 
total strain energy was less than 0.20 for the first 10 
modes (fig. 7(b)) and less than 0.28 for the succeeding 
10 modes. For the J-section frame, the strain-energy 
ratio in the web approached 0.4 in some of the modes 
(fig. 12(b)). 

7. For the I-section frame, the strain energy of 
the top flange is the dominant energy in the in- 
plane deformation modes, and the strain energy of 
the bottom flange (including the skin) dominates for 
the out-of-plane deformation modes. (See fig. 7(b).) 

8. The vibrational response of both the I-section 
and J-section frames is very sensitive to restraining 
the u\ displacements of the flanges (and skin). It is 
somewhat sensitive to the rotational restraint on <f> 2 . 
(See tables 4 and 5.) However, it is less sensitive to 
restraining the displacement components u 2 and w 9 
and the rotation (f >[ . 

9. The vibrational response of the I-section frame 
and J-section frame is more sensitive to variations 
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in the 45° or —45° fiber angles than to variations 
in the 0° or 90° fiber angles. The variations in the 
0° and 90° fibers of the web and the bottom flange 
have a noticeable effect on some of the modes, but 
their effect is generally less than that of the 45° and 
—45° fibers. (See figs. 9 and 14.) The vibrational 
response is also more sensitive to variations in the 
elastic moduli El and Glt than to any of the other 
material coefficients. (See figs. 10 and 15.) 

10. The sensitivity of the vibration frequencies 
with respect to variations in both El and Glt IS al- 
most the same for all the modes. (See figs. 10 and 15.) 
This uniform sensitivity may be attributed to the 
quasi-isotropic lamination used for both the flanges 
and the web. It suggests the feasibility of replacing 
the quasi-isotropic composite, in the one-dimensional 
thin-walled-beam model, with an equivalent isotropic 
material, as was done in the present study. 

Comments on Sources of Errors and 
Model Adjustment Techniques 

Sources of Errors 

The determination of natural frequencies and 
modes from vibration tests and numerical models in- 
volves numerous possible sources of discrepancies or 
errors that are related to mechanical and equipment 
limitations and to theoretical and physical assump- 
tions. The errors in vibration tests include inexact 
equipment calibration, excessive noise, manufactur- 
ing variations, incorrect transducer locations, and op- 
eration in a region of nonlinearity of the response. 
Numerical modeling errors can be attributed to in- 
accuracies in estimated material properties and to in- 
sufficient modeling detail. In the present study, care 
was exercised in collecting and recording the vibra- 
tion test data and in the selection of the numerical 
model. However, nominal material properties and 
lay-ups (fiber orientation of the different layers) were 
used in the numerical model. The sensitivity analysis 
helped identify the material and lamination parame- 
ters that need to be accurately determined. 

Model Adjustment Techniques 

In recent years, considerable efforts have been 
directed at improving and modifying the numerical 
model to obtain a better correlation with test results. 
These efforts started as trial-and-error approaches 
and evolved into systematic system identification 
and model adjustment techniques. Although these 
model adjustment techniques have not been used in 
the present study, the techniques are particularly 
useful for validating numerical models to be used in 
simulating transient dynamic response. 


Most of the model adjustment techniques are 
based on using the experimental modal data ( mea- 
sured eigenvalues and eigenvectors) to update the 
stiffness and/or mass matrices of the structure (e.g., 
see Berman 1979; Chen 1979; Wei 1980; Berman and 
Wei 1981; Baruch 1982; Grossman 1982; Berman 
and Nagy 1983; Jensen and Crawley 1984; Kabe 
1985; and Arruda and dos Santos 1989) and the two 
monographs (Ewins 1986 and Martinez and Miller 
1985). In some of the recent techniques, the sensitiv- 
ity derivatives with respect to the physical parame- 
ters of the numerical model are used in conjunction 
with optimization algorithms to obtain corrected (or 
adjusted) values of the physical parameters. 

Conclusions 

A detailed study is made of the effects of vari- 
ations in lamination and material parameters of 
thin-walled composite frames on their vibrational 
characteristics. The structures considered are semi- 
circular thin-walled frames with I- and J-section 
frames. The flanges, web, and skin of the stiffen- 
ers have quasi-isotropic laminations and the fiber 
orientation is made up of combinations of ±45°, 
0°, and 90° layers. Two computational models are 
used for predicting the vibrational characteristics. In 
the first model, the flanges and webs of the stiff- 
eners are modeled by using two-dimensional shell 
(and plate) finite elements. The second model is a 
finite-element discretization of the one-dimensional 
Vlasov-type thin-walled-beam theory. A mixed for- 
mulation is used with the fundamental unknowns 
consisting of both the generalized displacements and 
stress resultants (or internal forces) in the frame. 
The frequencies and modes predicted by the compu- 
tational models are compared with those obtained 
from experiments. A detailed study is made of 
the sensitivity of the vibration response to vari- 
ations in the fiber orientation, material proper- 
ties of the individual layers, and boundary condi- 
tions. On the basis of this study, the following 
conclusions are justified: 

1. For some of the higher vibration modes, the 
experimental frequencies for thin-walled frames are 
generally between those for the totally and partially 
clamped supports. 

2. Identification of the modes as in plane or out 
of plane can best be accomplished by examining the 
energy components associated with the in-plane and 
out-of-plane forces. Also, the minimum error to 
be expected when using one-dimensional thin-walled 
beams can be estimated by computing the ratio of 
the energy associated with the forces neglected in 
thin-walled beams to the total energy. 
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3. For quasi-isotropic composite frames, the vi- 
bration frequencies, associated with the lower modes, 
can be accurately predicted by an isotropic one- 
dimensional-beam model (with effective elastic mod- 
uli) . The accuracy of predictions is dependent on the 
cross-sectional distortions during the beam deforma- 
tions. As the cross-sectional distortions increase, the 
degradation of accuracy becomes more pronounced. 

4. The vibrational response of thin-walled semi- 
circular frames is very sensitive to restraining the 
u[ displacement component of the flanges along the 
length of the frame. It is somewhat sensitive to 
the restraint on the associated rotation component. 
However, it is less sensitive to restraining the other 
displacement and rotation components. 

5. The vibrational response of thin- walled com- 
posite frames with quasi- isotropic laminations is 
more sensitive to variations in the +45° or —45° fiber 
angles than to variations in the 0° or 90° fiber angles. 


Variations in the 0° and 90° fibers of the web and the 
bottom flange have a noticeable effect on some of the 
modes, but their effect is generally less than that of 
the 45° and —45° fibers. The vibrational response is 
also more sensitive to variations in the material co- 
efficients Ei and Git than to all other coefficients. 

6. The sensitivity of the vibration frequencies 
with respect to variations in the elastic moduli Ei 
and Git ls almost the same for all the modes be- 
cause of the quasi-isotropic lamination used for both 
the flanges and the web. It suggests the feasibility 
of replacing the quasi-isotropic composite with an 
equivalent isotropic material in the one-dimensional 
thin- walled-beam analysis, as was done in the present 
study. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
September 4, 1990 
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Appendix 

Fundamental Equations of 
Thin- Walled-Beam Theory Used in 
Present Study 

The fundamental equations of the linear, Vlasov- 
type theory of curved thin- walled beams are given in 
this appendix. A right-handed orthogonal coordinate 
system is used with the x-axis passing through the 
centroids of the cross sections. (See fig. 3.) The beam 
is assumed to be curved in one direction only (in the 
xz- plane). 

Displacement Assumptions 

Based on the assumption that the projection of 
each cross section on a plane normal to the initial 
centroidal axes does not distort during deformation, 
the displacement field in the plane of the cross section 
(yz- plane) is represented by 


z-directions, is the twist, and 7^ and 7^ are 
the transverse shear strains. The strain parameters 
7 xz> 7 Sir «* and axe functions 
of x only and can be expressed in terms of the 
displacement and rotation components as follows: 


0 q jlO , Vx 
*y = VVz + 

K 0 Z = 8(Py 
7 xy = dv° - 4>°z 

7 ° xz = ~+dw° + 4> 0 y 

jfi 

^0 _ , <Pz 

K t ~ d( Px + ^ 

= d0° 


(A3) 



where d = d/dx and R is the radius of curvature 
of the centerline of the beam. Also, the following 
constraint condition is used to relate 0° and (f)®: 

d<t>l -e° = o (A4) 


Constitutive Relations 

The relations betwen the internal forces and the 
strain components are given by 


where and are the axial and transverse 

displacement components at y = z = 0; (f)®, <f>y, and 
<)P Z are the rotation components about the coordinate 
axes; 0° is the rate of twist of the beam; and TD is 
the sectorial coordinate (warping of the cross section 
for a unit rate of twist). The seven generalized 
displacement parameters u°, t>°, iu°, <jP z , and 

0° are functions of x only. 

Strain Assumptions 

The following expressions are used for the three 
nonzero components of the strain field in the plane 
of the cross section: 


e x = £x~ y*y + 2*2 - ' 

7 xy = 7 xy ~ ZK t 
Ixz = 7 xz + l /«t 


(A2) 


where e % is the extensional strain of the centerline, 
K,y and k® are the curvature changes in the y- and 



10 



where A is the cross-sectional area; I y , I z , and I yz 
are the second moments of the cross section (mo- 
ments and product of inertia); J is the Saint- Venant 
torsion constant; I u is the principal second secto- 
rial moment of the cross section (sectorial moment 
of inertia); E and G are the effective Young’s and 
shear moduli of the material; N x is the axial force; 
M y and M z are the bending moments; B u is the 
bimoment; Q y and Q z are the transverse shearing 
forces; and Mt is the twisting moment. The defini- 
tion of the sectorial properties of the cross section is 
given in Vlasov (1961), Zbirohowski-Koscia (1967), 
and Gjelsvik (1981). 


Variational Functional 

The functional used in the element development 
is given by 


* = *HR + 




L 


i „ „ 

(A) 2 dx (A7) 


where A is the Lagrange multiplier, e is a penalty 
parameter, l is the length of the element, and khr 
is the functional of the Hellinger-Reissner mixed 
variational principle. The expression for ttjjr is 


*HR= f\v-U c + K) dx 
Jo 


(A8) 


where 


V 


( e 


N x 

M z 

My 

Bfjj / 


> < 


0 X 


l J 


> + 


Qy 

Qz 


' Ixy 1 

< 7 2, 


(A9) 



■ N x ' 

t 

■A ■ 

1 

M z 

> 

Iz ~Iyz 

2 E < 

My I 


Iy ■ 


Ib u J 


Symm I w . 


-1 

(Nx) 


M z 

< 

\M y \ 

- 

l B u ) 





L Symm 


7J 



(A10) 


K = | w 2 (a [(u 0 ) 2 + ( v °) 2 + (w 0 ) 2 ] 

+ {Iy + Iz){4>x) 2 + lyi^ y) 2 
+ I z (<p° z ) 2 - 2 I yz <t>y z + /„(0 0 ) 2 } (All) 


where p is the mass density of the material. In 
equations (A5) to (All), y and 2 are centroidal 
coordinates. (See fig. 3.) A Fortran program for 
evaluating the principal sectorial properties is listed 
by Coyette (1987). 
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Table 1 . Measured Thicknesses 


6, deg 

/tl, in. 

/i 2 > in - 

h 3 , in. 

/14 j in . 

/ 15 , in. 

/i 6 , in. 

2 

0.046 

0.042 

0.051 

0.133 

0.139 

0.089 

10 

.041 

.042 

.046 

.119 

.125 

.078 

20 

.042 

.043 

.039 

.127 

.130 

.084 

30 

.041 

.043 

.041 

.127 

.126 

.085 

40 

.043 

.044 

.041 

.121 

.119 

.080 

50 

.045 

.046 

.046 

.120 

.116 

.081 

60 

.044 

.045 

.049 

.124 

.127 

.085 

70 

.043 

.044 

.055 

.124 

.131 

.080 

80 

.045 

.044 

.056 

.125 

.122 

.085 

90 

.051 

.049 

.050 

.128 

.129 

.088 

100 

.046 

.045 

.040 

.119 

.118 

.083 

110 

.045 

.047 

.039 

.121 

.122 

.081 

120 

.044 

.047 

.040 

.118 

.120 

.083 

130 

.045 

.045 

.041 

.118 

.122 

.083 

140 

.045 

.044 

.042 

.120 

.124 

.082 

150 

.044 

.043 

.041 

.120 

.120 

.085 

160 

.043 

.043 

.040 

.123 

.125 

.086 

170 

.043 

.042 

.052 

.120 

.117 

.085 

178 

.042 

.044 

.048 

.136 

.141 

.093 


Dimensions for I-Section Frame 


/ 17 , in. 

h, in- 

62 , in. 

63, in- 

64 , in. 

di, in. 

d, 2 , in. 

0.092 

0.712 

0.635 

1.790 

1.625 

0.785 

0.770 

.086 

.800 

.633 

1.795 

1.670 

.785 

.785 

.092 

.800 

.653 

1.775 

1.670 

.800 

.800 

.091 

.773 

.680 

1.780 

1.680 

.800 

.810 

.082 

.775 

.685 

1.780 

1.687 

.815 

.810 

.081 

.780 

.675 

1.770 

1.685 

.820 

.810 

.088 

.795 

.685 

1.790 

1.690 

.820 

.815 

.093 

.825 

.680 

1.775 

1.680 

.825 

.810 

.086 

.825 

.650 

1.767 

1.687 

.815 

.810 

.091 

.850 

.645 

1.780 

1.660 

.828 

.790 

.080 

.830 

.645 

1.800 

1.675 

.812 

.800 

.081 

.840 

.635 

1.790 

1.690 

.810 

.790 

.079 

.800 

.670 

1.770 

1.687 

.810 

.780 

.078 

.755 

.690 

1.765 

1.705 

.810 

.790 

.079 

.730 

.725 

1.770 

1.675 

.805 

.795 

.079 

.730 

.708 

1.785 

1.670 

.805 

.785 

.082 

.710 

.712 

1.810 

1.655 

.800 

.808 

.077 

.700 

.728 

1.795 

1.655 

.807 

.785 

.096 

.692 

.700 

1.810 

1.650 

.800 

.785 



Table 2. Measured Thicknesses and Dimensions for J-Section Frame 




















Table 3. Decomposition of Total Complementary Strain Energy Into Components 


[U c = u l + u 2 + u 3 \ 



Associated stress resultants 
(see fig. 2) 


Energy 

components 

Web 

Flanges and skin 

Comments 

U i 

N u N U 

Ni, M\,Q\ 

In-plane response quantities 

U 2 

Mi, M 12 , <?i 

N\ 2 , Mi2 

Out-of-plane response quantities 

Uz 

n 2 , m 2 , q 2 

n 2 , m 2 , q 2 

Response quantities neglected in 
one-dimensional model 




Table 4. Effect of Boundary Conditions on Frequencies Obtained by Two-Dimensional 
Finite-Element Model for I-Section Frame 


Numbers in parentheses refer to ratios of partially clamped 
to totally clamped model frequencies 


Mode 

Frequencies 
of totally 
clamped 
model, 
Hz 

Frequencies of partially clamp 
(with the following genera 
displacements unrestrainec 

ed model 
lized 
0) Hz 

Experimental 

frequencies, 

Hz 

<t>2 

<t>2 and <f>3 

u[ in 
flanges 

u[ in flanges 
and <f > 2 

1 

9.201 

9.001 

9.001 

6.788 

6.632 

9.2 



(0.978) 

(0.978) 

(0.738) 

(0.721) 


2 

31.86 

31.06 

31.06 

18.11 

17.87 

29.7 



(0.975) 

(0.975) 

(0.568) 

(0.561) 


3 

37.52 

37.37 

37.37 

34.17 

33.36 

35.9 



(0.996) 

(0.996) 

(0.911) 

(0.889) 


4 

73.85 

71.82 

71.81 

38.09 

37.69 

66.6 



(0.973) 

(0.972) 

(0.516) 

(0.510) 


5 

81.34 

81.03 

81.03 

74.30 

73.44 

78.1 



(0.996) 

(0.996) 

(0.913) 

(0.903) 


6 

133.9 

130.1 

130.1 

75.56 

74.14 

119.0 



(0.972) 

(0.972) 

(0.564) 

(0.554) 


7 

149.2 

148.6 

148.6 

129.6 

128.1 

145.0 



(0.996) 

(0.996) 

(0.869) 

(0.858) 


8 

203.3 

198.1 

198.1 

139.8 

137.4 

193.0 



(0.974) 

(0.974) 

(0.688) 

(0.676) 


9 

226.5 

225.6 

225.6 

199.3 

196.9 

216.0 



(0.996) 

(0.996) 

(0.880) 

(0.870) 

223.0 

10 

281.9 

275.2 

275.2 

214.2 

210.9 

260.0 



(0.976) 

(0.976) 

(0.760) 

(0.748) 


11 

320.6 

319.3 

319.3 

268.0 

264.9 

309.0 



(0.996) 

(0.996) 

(0.836) 

(0.826) 


12 

349.8 

342.3 

342.3 

305.1 

300.7 

(Missed) 



(0.979) 

(0.979) 

(0.872) 

(0.860) 


13 

419.1 

412.6 

412.6 

343.0 

339.0 

401.0 



(0.985) 

(0.985) 

(0.819) 

(0.809) 
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Table 5. Effect of Boundary Conditions on Frequencies Obtained by Two-Dimensional 
Finite-Element Model for J-Section Frame 

Numbers in parentheses refer to ratios of partially clamped 
to totally clamped model frequencies 


Mode 

Frequencies 
of totally 
clamped 
model, 
Hz 

Frequencies of partially clamped model 
(with the following generalized 
displacements unrestrained), Hz 

Experimental 

frequencies, 

Hz 

<t>2 

<t >2 and <t>' 3 

u[ in 
flanges 

u[ in flanges 
and <f > 2 

i 

11.53 

11.24 

11.24 

8.488 

8.408 

11.6 



(0.975) 

(0.975) 

(0.736) 

(0.729) 


2 

36.87 

36.64 

36.64 

22.41 

22.37 

32.1 



(0.994) 

(0.994) 

(0.608) 

(0.607) 


3 

39.81 

38.80 

38.79 

32.77 

32.18 

37.0 



(0.975) 

(0.974) 

(0.823) 

(0.808) 


4 

79.22 

78.99 

78.99 

48.32 

47.94 

69.0 



(0.997) 

(0.997) 

(0.610) 

(0.605) 


5 

91.41 

88.81 

88.78 

72.64 

71.68 

79.0 



(0.972) 

(0.971) 

(0.795) 

(0.784) 


6 

143.9 

143.5 

143.5 

96.58 

95.35 

126.0 



(0.997) 

(0.997) 

(0.671) 

(0.663) 


7 

168.1 

163.6 

163.5 

134.4 

132.9 

145.0 



(0.973) 

(0.973) 

(0.800) 

(0.790) 


8 

214.1 

213.4 

213.4 

167.2 

165.0 

191.0 



(0.997) 

(0.997) 

(0.781) 

(0.777) 


9 

263.0 

256.9 

256.8 

206.5 

204.3 

221.0 



(0.977) 

(0.976) 

(0.785) 

(0.770) 

229.0 







247.0 

10 

297.6 

296.2 

296.2 

251.5 

248.3 

266.0 



(0.995) 

(0.995) 

(0.845) 

(0.834) 


11 

368.2 

361.2 

361.1 

298.7 

295.5 

339.0 



(0.981) 

(0.981) 

(0.811) 

(0.803) 


12 

382.8 

380.3 

380.3 

336.7 

332.8 

347.0 



(0.993) 

(0.993) 

(0.880) 

(0.869) 


13 

468.2 

462.0 

461.9 

402.4 

398.0 

403.0 



(0.987) 

(0.987) 

(0.859) 

(0.850) 










Flanges and skin 



Figure 2. Sign convention for generalized displacements and stress resultants in two-dimensional model. 



Figure 3. Sign convention for generalized displacements and stress resultants in one-dimensional model. 
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L-89-4693 

Figure 4. Thin- walled semicircular graphite-epoxy specimens and equipment. 
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Frame section 



Figure 5. Schematic of test apparatus. 
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(a) Results for experimental and two-dimensional model. 



(b) Results for experimental and bounding two-dimensional model. 

Figure 6. Comparison of finite-element and experimental frequencies for thin-walled composite frame with 
I cross section. 
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Figure 6. Concluded. 
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Figure 7. 


Energy components in different vibration modes of thin-walled composite frame with I cross 


section. 
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26 


Figure 8. Mode shapes associated with lowest ten frequencies for thin-walled composite beam with I cross section. Numbers in parentheses 
are experimental frequencies. Other numbers are frequencies obtained by clamped finite-element model. The plots in the last two columns 
are not drawn to scale. 
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Figure 8. Concluded. 



(a) Top flange. 



(b) Web. 

Figure 9. Sensitivity of vibration frequencies to variations and fiber orientation of flanges and web for 
thin-walled composite frame with I cross section. 
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Mode 

(c) Bottom flange and skin. 
Figure 9. Concluded. 
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Figure 10. Sensitivity of vibration frequencies to variations in material characteristics of thin-walled 
composite frames with I cross section. 
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(a) Results for experimental and two-dimensional model. 



(b) Results for experimental and bounding two-dimensional model. 

Figure 11. Comparison of finite-element and experimental frequencies for thin- walled composite frame 
with J cross section. 
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Figure 11. Concluded. 
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Mode 


(a) Ui, U 2 , and U 3 . 



Mode 

(b) fAu? and 

Figure 12. Energy components in different vibration modes of thin-walled composite frame with J cross 
section. 
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Figure 13. Mode shapes associated with lowest ten frequencies for thin-walled composite beam with J cross section. Numbers in parentheses 
are experimental frequencies. Other numbers are frequencies obtained by clamped finite-element model. The plots in the last two columns 
are not drawn to scale. 
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Figure 13. Concluded. 





(a) Top flange. 
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(b) Web. 


Figure 14. Sensitivity of vibration frequencies to variations and fiber orientation of flanges and web for 
thin-walled composite frame with J cross section. 
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(c) Bottom flange and skin. 
Figure 14. Concluded. 
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Figure 15. Sensitivity of vibration frequencies to variations in material characteristics of thin-walled 
composite frames with J cross section. 
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